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Supersonic beams of polar molecules are deflected using inhomogeneous electric fields. The quantum- 
state selectivity of the deflection is used to spatially separate molecules according to their quantum 
state. A detailed analysis of the deflection and the obtained quantum-state selection is presented. 
The rotational temperatures of the molecular beams are determined from the spatial beam profiles 
and are all approximately 1 K. Unprecedented degrees of laser-induced alignment ((cos^^2d) = 
0.972) and orientation of iodobenzene molecules are demonstrated when the state-selected samples 
are used. Such state-selected and oriented molecules provide unique possibilities for many novel 
experiments in chemistry and physics. 

PACS numbers: 37.20.+j, 33.15.-e 



I. INTRODUCTION 

For a large range of experiments in chemistry and 
physics, a high level of control over the external and 
internal degrees of freedom of molecules is very benefi- 
cial. This includes control over the translational and the 
rotational motions, as well as the selection of a single 
quantum state or a small set of states. Such quantum- 
state-selected targets provide unique possibilities, for ex- 
ample, for manipulating the external degrees of freedom 
with static electric fields [1, 2] or optical fields [3, 4], 
or both [5, 6]. The quantum- state selection also natu- 
rally discriminates between individual stereo-isomers of 
large molecules [7]. The resulting samples of aligned 
or oriented individual isomers offer unique prospects for 
novel experiments with complex molecules, such as fem- 
tosecond pump-probe measurements, x-ray or electron 
diffraction in the gas-phase [8, 9], high-harmonic gener- 
ation [10], or tomographic reconstructions of molecular 
orbitals [11]. Moreover, it would provide considerably 
increased control in reaction dynamics experiments [12]. 

Strong cooling can be achieved in supersonic expan- 
sions of molecules seeded in an inert atomic carrier gas. 
For small molecules (consisting of just a few atoms) only 
a few rotational states are populated at the typical tem- 
peratures on the order of 1 K. For larger polyatomic 
systems rotational cooling down to or even below 1 K 
still leaves the molecular ensemble distributed over a 
considerable number of rotational states, thereby often 
masking quantum-state- specific effects. State selection 
can be performed using inhomogeneous electric or mag- 
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netic fields. The possibility to deflect polar molecules 
in a molecular beam with an electric field was first de- 
scribed by Kallmann and Reiche in 1921 [13] and exper- 
imentally demonstrated by Wrede in 1927 [14]. As early 
as 1926, Stern suggested that the technique could be 
used for the quantum-state separation of small diatomic 
molecules at low temperatures [15]. In 1939 Rabi intro- 
duced the molecular beam resonance method, by using 
two deflection elements of oppositely directed gradients 
in succession, to study the quantum structure of atoms 
and molecules [16]. Whereas deflection experiments al- 
low the spatial dispersion of quantum states, they do not 
provide any focusing. For small molecules in low-field- 
seeking states this issue could be resolved using multi- 
pole focusers with static electric fields. These were de- 
veloped independently in 1954/55 in Bonn [17, 18] and in 
New York, where they were used to produce the popula- 
tion inversion for the first MASER experiments [19, 20]. 
About ten years later, molecular samples in a single ro- 
tational state were used for state specific inelastic scat- 
tering experiments by the Bonn group [21] and, shortly 
thereafter, for reactive scattering [22, 23]. In the follow- 
ing decades, multipole focusers were extensively used to 
study steric effects in gas-phase reactive scattering exper- 
iments [12, 24]. Also for the investigation of steric effects 
in gas-surface scattering [25] and photodissociation [26] 
experiments, the preparation of oriented samples of state- 
selected molecules using electrostatic focusers was essen- 
tial. For about ten years, it is also possible to manip- 
ulate the speed of small molecules using switched inho- 
mogeneous electric fields in the so-called Stark deceler- 
ator [27]. More recently, also its optical [28] and mag- 
netic [29] analogs have been demonstrated. 

Obtaining similar control over large molecules is more 
difficult, because all low-lying quantum states are high- 
field seeking at the required electric field strengths. 
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In order to confine these molecules, dynamic focusing 
schemes are necessary [30-32] . Dynamic focusing of large 
molecules has been demonstrated in the alternating- 
gradient (AG) deceleration of benzonitrile [33] and in the 
conformer selection of 3-aminophenol [7]. However, if fo- 
cusing is not necessary, spatial dispersion of quantum 
states can still be achieved using static electric fields in 
a Stern- Gerlach- type deflector. This molecular beam de- 
flection has been used extensively as a tool to determine 
dipole moments and polarizabilities of molecular systems 
ranging from diatomics [14] over clusters [34] to large 
biomolecules [35]. 

Recently, we have demonstrated the quantum-state se- 
lection of large molecules [36] following the original pro- 
posal of Stern [15]. Here, the details of the electrostatic 
deflection are presented. It is shown, how the rotational 
temperature of cold supersonic jets can be determined 
with high precision from deflection measurements and 
that indeed a small subset of quantum states can be ad- 
dressed in deflected samples of large molecules. In par- 
ticular, the ground state has the largest Stark shift and 
molecules residing in this state are deflected the most. 
Our goal is to isolate and use rotational ground state 
molecules, or at least samples of molecules in the few 
lowest lying states, as targets for various experiments. 
Since the deflection does not change the initial state dis- 
tribution but merely disperses it, it is crucial that the 
population of ground state molecules in the molecular 
beam is initially as large as possible. Therefore, the ro- 
tational temperature of the molecular beam is made as 
low as possible using a high-pressure supersonic expan- 
sion [37]. It is shown how the resulting state-selected 
molecules can be used to improve one-dimensional (ID) 
laser-induced alignment [4, 38] and mixed field orienta- 
tion [36, 39, 40]. Here, alignment refers to confinement of 
a molecule-fixed axis (typically, the largest polarizability 
axis) along a laboratory-fixed axis and orientation refers 
to the molecular dipole moments pointing in a partic- 
ular direction. Alignment and orientation occur in the 
adiabatic limit where the laser field, used to align the 
molecules, is turned on and off slowly compared to the 
inherent rotational periods of the molecule [41, 42]. The 
state selection leads to strong enhancement in the degree 
of orientation and alignment of iodobenzene molecules 
compared to that achieved when no deflection is used. 



II. EXPERIMENTAL 

A schematic of the experimental setup is shown in Fig- 
ure 1. The molecular beam machine consists of three dif- 
ferentially pumped vacuum chambers; The source cham- 
ber housing a pulsed valve (pumped by a 2000 1/s tur- 
bomolecular pump), the deflector chamber (pumped by 
a 500 1/s turbomolecular pump) and the detection cham- 
ber housing the ion/electron spectrometer (pumped by a 
500 1/s turbomolecular pump). About 3 mbar of iodoben- 
zene (Sigma Aldrich, 98 % purity) or benzonitrile (Sigma 
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FIG. 1: (Color online): Scheme of the experimental setup. In 
the inset, a cut through the deflector is shown, and a contour- 
plot of the electric field strength is given. Details of the ve- 
locity map imaging spectrometer are shown in Figure 7. See 
text for details. 



Aldrich, 98 % purity) is seeded in an inert carrier gas 
and expanded through a pulsed valve into vacuum. In 
order to obtain optimal cooling of the molecular beam, 
a miniaturized, high pressure Even-Lavie valve [37] is 
used operating at a backing pressure of 90 bar of He 
or 20 bar of Ne, limited by the onset of cluster forma- 
tion. While rotational temperatures down to 0.4 K have 
been achieved under similar conditions [43], the typi- 
cal rotational temperature in our experiments is ~ 1 K. 
Two 1-mm-diameter skimmers placed 15 cm (separating 
the source and the deflector chamber) and 38 cm down- 
stream from the nozzle collimate the molecular beam be- 
fore it enters a 15-cm-long electrostatic deflector. A cut 
through the electrodes of the deflector is shown in the 
inset of Figure 1 together with the electric field created. 
A trough with an inner radius of curvature of 3.2 mm 
at ground potential and a rod with a radius of 3.0 mm 
at high voltage create a two- wire field [44]. The verti- 
cal gap across the molecular beam axis is 1.4 mm, while 
the smallest distance between the electrodes is 0.9 mm. 
The two-wire field geometry is ideally suited for molec- 
ular beam deflection. The gradient of the electric field 
along the vertical direction is large and nearly constant 
over a large area explored by the molecular beam, while 
the electric field is very homogeneous along the horizon- 
tal direction. Thus, a polar molecule experiences a nearly 
constant force in the vertical direction independent of its 
position within the deflector, while the force in the hori- 
zontal direction (i. e., broadening of the beam in the hor- 
izontal direction) is minimized. In our setup, the deflec- 
tor is mounted such that molecules in high-field-seeking 
(low- field-seeking) quantum states are deflected upwards 
(downwards). 

After passing through the deflector, the molecular 
beam enters the differentially pumped detection chamber 
via a third skimmer of 1.5 mm diameter. In the detection 
area, the molecular beam is crossed by one or two laser 
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beams that are focused by a spherical lens with a focal 
length of / = 300 mm. The lens is mounted on a verti- 
cal translation stage so that the height of the laser foci 
can be adjusted with high precision. In the first part of 
the experiment, where the beam deflection of iodoben- 
zene and benzonitrile is characterized, only one laser, 
the probe laser, is used. This Ti: Sapphire laser (25 fs 
(FWHM) pulses, 800 nm, beam-waist coq = 21 pm) is 
used to determine the relative density in the molecular 
beam via photoionization. In the second part of the ex- 
periment, an additional laser pulse is included to study 
laser-induced alignment and orientation of iodobenzene. 
For these experiments, 10 ns (FWHM) long pulses from 
a NdiYAG laser (1064 nm, ujq = 36 pm) are overlapped 
in time and space with the probe laser pulses. While 
the YAG laser induces adiabatic alignment and orienta- 
tion, here the fs-laser is used to determine the spatial 
orientation of the target molecules via Coulomb explo- 
sion. Ionic fragments produced in the Coulomb explosion 
are accelerated in a velocity focusing geometry towards 
the detector. The detector can be gated with a time 
resolution of ~ 90 ns, which allows for mass selective 
detection of individual fragments. A microchannel plate 
(MCP) detector backed by a phosphor screen is employed 
to detect the position of mass-selected ions. In particu- 
lar, 1+ fragment ions, formed in the Coulomb explosion 
of iodobenzene, are particularly useful experimental ob- 
servables since they recoil along the C-I symmetry axis 
of the molecule. Thus, 2D ion images of 1+ recorded with 
a CCD camera provide direct information about the in- 
stantaneous molecular orientation of the C-I bond axis 
with respect to the laboratory frame and are, therefore, 
the basic observables in these experiments. All experi- 
ments are conducted at 20 Hz, limited by the repetition 
rate of the YAG laser. 



III. RESULTS AND DISCUSSION 

In the first experiments shown in Section HI A, a de- 
tailed analysis of electrostatic beam deflection is pre- 
sented. It is shown, that electrostatic deflection of cold 
molecular beams can be used to determine the rotational 
temperature of a supersonic jet. Furthermore, the de- 
gree of quantum-state selectivity that can be achieved 
with our setup is investigated. As an application it is 
demonstrated, how this quantum-state selectivity of the 
deflection process can be exploited to obtain an unprece- 
dented degree of laser-induced alignment (Section IIIB) 
and orientation (Section HI C) of iodobenzene molecules. 



A. Electrostatic Deflection of Cold Molecular 
Beams 

In the first experiment, the deflection of benzonitrile 
molecules (BN, C7H5N) seeded in 90 bar of He is in- 
vestigated. BN is an ideal candidate for electrostatic 
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FIG. 2: (Color online): Energy as a function of the elec- 
tric field strength for the lowest rotational quantum states 
of benzonitrile. The shaded area indicates the actual electric 
field strength inside the deflector. In the inset, the molecu- 
lar structure is shown together with the relevant molecular 
constants [45]. 

beam deflection due to its large permanent dipole mo- 
ment of 4.515 D. From the precisely known molecular 
constants [45] the energy of a given rotational quantum 
state can be calculated as a function of the electric field 
strength. The exact procedure is detailed in Section A 1. 
Figure 2 shows the Stark energies for the lowest rota- 
tional states of BN. Due to the small rotational constants 
and the resulting high density of rotational states, a large 
number of states is populated even under the cold condi- 
tions in a supersonic expansion. At a rotational temper- 
ature of 1 K, the typical temperature in our experiments 
{vide infra) ^ 66 rotational quantum states (with 419 M- 
components) have a population larger than 1 % relative 
to the ground state. At the electric field strengths present 
in the deflector, indicated by the shaded area in Figure 2, 
all low-lying quantum states are high-field seeking. This 
is due to mixing of closely spaced states of the same sym- 
metry and is typical for large asymmetric top molecules. 
The Stark shift and thus the force a molecule experiences 
in an inhomogeneous electric field depends on the rota- 
tional quantum state. Molecules in the ground state have 
the largest Stark shift and are, therefore, deflected the 
most. In general, the Stark shift decreases with increas- 
ing J quantum number. Thus, the lower the rotational 
temperature of the molecular beam is, the more the beam 
is deflected. 

Figure 3 shows vertical intensity profiles of BN for var- 
ious high voltages applied to the deflector. Vertical in- 
tensity profiles are obtained by recording the BN+ signal 
from photoionization by the fs-laser as a function of the 
vertical position of the laser focus. If no high-voltage 
is applied to the deflector, the molecular beam extends 
over about ^ 2 mm. In this case, the size of the molecular 
beam in the detection region is determined by the me- 
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FIG. 3: (Color online): The vertical spatial profile of the 
molecular beam for different deflection voltages applied, mea- 
sured by recording the laser-induced BN^ signal (see text). 
The experimental data are shown as symbols together with 
the corresponding simulated profiles (lines). 



FIG. 4: (Color online): The vertical profile of the molecular 
beam measured by recording the laser-induced BN^ signal 
(see text). The experimental data are shown together with 
the corresponding simulated profiles (solid lines) for different 
rotational temperatures. Simulations are shown for the two- 
temperature model; see text for details. 



chanical aperture of the experimental setup, i.e., by the 
dimensions of the deflector and the last skimmer before 
the detection region. As the high voltage is turned on, 
the molecular beam profile broadens and shifts upwards. 
At a voltage of 10 kV, a large fraction of the molecules 
is deflected out of the original, undeflected beam profile. 
A small fraction of the molecules, however, is almost un- 
affected by the deflector. 

In order to understand these experimental findings, 
Monte Carlo simulations are employed, which are de- 
scribed in detail in Section A 2. In brief, trajectory 
calculations are performed for molecular packets of in- 
dividual rotational quantum states. These calculations 
yield single-quantum-state deflection profiles. Then, the 
single-state profiles are averaged according to the pop- 
ulations of the respective states in the original molecu- 
lar beam (i.e., at the entrance of the deflector). From 
these simulations it is obvious, that the molecules in 
the original beam are not rotationally thermalized, an 
effect that has previously been observed in rotation- 
ally resolved spectroscopy [46, 47]. A number of differ- 
ent descriptions of the populations of rotational states 
have been given [47]; we use the formula for a two- 
temperature model originally proposed by Levy and 
coworkers [46]. For details and the approximation of the 
high-temperature component see Section A3. Finally, 
the rotational temperature of the low-temperature com- 
ponent in the molecular beam is obtained by fitting the 
simulated deflection profiles to the experimental data us- 
ing a local optimization algorithm. All deflection pro- 
files measured at the different voltages are fitted simul- 
taneously, where the fraction q of the low-temperature 
component, a general intensity scaling factor s of the 
deflected profiles (with respect to the undeflected beam 



profile), and the rotational temperature Trot of the low- 
temperature component are the fitting parameters. Best 
agreement between experimental data and simulations is 
found for = 0.93 and Trot = 0.8 K. The resulting simu- 
lated deflection profiles nicely reproduce the experimen- 
tal data as shown in Figure 3 (solid lines). In particular, 
the undeflected part of the molecular beam for 10 kV 
can be perfectly simulated, which indicates that the use 
of a two-temperature model was indeed justified. For 
comparison, also a simulated deflection profile for 10 kV 
using a one-temperature model is shown (dashed line in 
Figure 3). 

In order to estimate the uncertainty of Trot, deflection 
profiles are calculated for different rotational tempera- 
tures. For each fixed rotational temperature, the best 
values for s and q are determined using the fitting pro- 
cedure outlined above and the resulting deflection profile 
for a voltage of 10 kV is plotted in Figure 4. With in- 
creasing Trot, the peak of the beam profile shifts towards 
smaller y- values, while, at the same time, the intensity in 
the undeflected part of the beam profile is reduced. From 
the comparison of experimental data and simulation, an 
uncertainty of Trot of ±0.2 K is estimated. 

The deflection of iodobenzene molecules (IB, CeHsI) is 
investigated in the same way. Vertical intensity profiles 
are measured by recording the signal of 1+ ions, created 
by Coulomb explosion with a circularly polarized probe 
pulse, as a function of the vertical position of the probe 
laser focus. Figure 5 a) shows deflection measurements 
for IB seeded in 90 bar of He. IB (mass 204 u) is heavier 
than BN (mass 103 u) and has a considerably smaller 
dipole moment of only 1.625 D [48] compared to BN. 
Therefore, smaller deflection amplitudes are observed for 
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FIG. 5: (Color online): The vertical profile of the molecu- 
lar beam for different deflection voltages applied, measured 
by recording the laser-induced signal (see text). The ex- 
perimental data are shown together with the corresponding 
simulated profiles. Figure (a) shows the beam profiles for 
iodobenzene seeded in 90 bar He, and Figure (b) shows the 
respective profiles for iodobenzene seeded in 20 bar Ne. Ar- 
rows indicate the laser positions for alignment and orientation 
experiments {vide infra). 



IB under identical expansion conditions. However, the 
interaction time with the electric field and the time-of- 
fiight from deflector to detection region can be increased 
when Ne is used as a carrier gas instead of He. Changing 
the carrier gas reduces the mean velocity of the molec- 
ular beam from ~ 1800 m/s to ~ 800 m/s and signifi- 
cantly enhances the observed deflection as shown in Fig- 
ure 5 b). Following the fitting procedure outlined above, 
Trot can be determined for IB as well. In the case of 
IB, a one-temperature model with a rotational temper- 
ature of 1.05 K fits the experimental data best for the 
deflection measurements in helium as well as in neon. 
The uncertainty of Trot is estimated to be ±0.1 K for IB 
seeded in Ne and ±0.2 K for IB in He. The somewhat 
larger uncertainty for the measurements in He reflects 
the small deflections observed for He. The simulated de- 
flection profiles for IB are shown as solid lines in Figure 5 
and agree well with the experimental data. 

The main purpose of the deflection studies presented 
in this work is to provide quantum-state-selected samples 
of large molecules for further experiments. The degree of 
deflection that a molecule experiences in the electric field 
of the deflector depends on its quantum state. The rel- 
evant quantity is the effective dipole moment /ieff (the 
negative slope of the Stark curve), which depends on the 
electric field strength. Molecules residing in low rota- 
tional quantum states have generally the largest /ieff and 
are, therefore, deflected most. These molecules can sim- 
ply be addressed by moving the laser focus in the de- 
tection region towards the upper cut-off of the molecular 
beam profile. In order to understand the laser-induced 
alignment and orientation experiments presented in Sec- 
tion HI B and HI C, it is crucial to know the relative pop- 
ulations of individual quantum states that are probed at 
a given height of the laser focus. The positions of the 
laser foci within the molecular beam profile during the 
alignment and orientation measurements are indicated 
by arrows in Figure 5. At this position, the intensity 
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FIG. 6: (Color online): (a) Energy as a function of the electric 
field strength for selected quantum states of IB. Solid lines 
represent quantum states that are present close to the upper 
cut-off of the molecular beam profile for IB seeded in Ne. The 
respective JxaK^ quantum numbers are given in the figure, 
(b) Effective dipole moment for selected quantum states of IB. 
The shaded area represents the range of electric field strengths 
in the deflector at 10 kV. 



of the deflected beam is 9 % of the undeflected peak 
intensity. The composition of the molecular packets in 
the detection region can be extracted from the simulated 
deflection profiles. Table I provides an overview of the 
most abundant quantum states present in different re- 
gions of the beam profiles for IB at a rotational temper- 
ature of 1.05 K. For comparison, also the population of 
each rotational quantum state in the undeflected beam 
for this rotational temperature is given. At the posi- 
tion of the laser focus for the orientation experiments 
in Ne (column four in Table I), the population of the 
lowest quantum states is significantly enhanced in the 
deflected beam compared to the undeflected beam. The 
fraction of ground state molecules is enhanced by a fac- 
tor of five, for instance. About 97 % of the population 
resides in the quantum states listed in Table I with the 
jRaKc = 3o3 state being most abundant. Moving the 
laser focus even closer towards the upper cut-off in the 
beam profile should reduce the number of quantum states 
that are probed even further. If a reduction of the beam 
intensity by two orders of magnitude (compared to the 
undeflected beam) can be afforded, only 4 quantum states 
are predicted to be probed with 37 % of the molecules 
being in the JxaK^ = ^03 state. At first glance, it is sur- 
prising that this state and not the absolute ground state, 
which is expected to have the largest /ieff 7 is populated 
most in the deflected beam. In order to understand this, 
the Stark curves for the most deflected quantum states 
of IB are shown in Figure 6 a), together with their ef- 
fective dipole moments (Figure 6 b). Below the relevant 
electric field strengths, both the JxaK^^ = 3o32 and 
the JRaKc^ = 4i33 M-sublevels have avoided crossings 
with close-by states of the same symmetry (dashed lines 
in Figure 6 a). These avoided crossings lead to large lo- 
cal effective dipole moments that are comparable to the 
ground- state /ieff- [79] Thus molecules in these quantum 
states are deflected as much as ground-state molecules. 
Furthermore, states with M 7^ are doubly degenerate, 
whereas the ground state with M = is only singly de- 
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TABLE I: Relative population of individual quantum states in the deflected part of the molecular beam profile for Trot = 1-05 K. 
Left: IB in He at 1 % of peak intensity of undeflected beam. Center: IB in Ne at 1 % of peak intensity of undeflected beam. 
Right: IB in Ne at 9 % of peak intensity of undeflected beam (here orientation images were taken). Pm denotes relative 
population of individual M-sublevels in %, Pj^^^j^^ the sum over all M-sublevels, and Pj^^j^ the relative population of a given 
rotational quantum state in a free jet. 



generate. Therefore, the population of molecules in the 
jRaKc^ — 3o32 in the undeflected beam is already larger 
than the population in the ground state. [80] From Ta- 
ble I it is clear, that it will be difficult to isolate the rota- 
tional ground state of iodobenzene in our setup. Never- 
theless, given that the fraction of ground-state molecules 
could be increased from 1 % in the undeflected to 26 % in 
the deflected beam, dramatic effects are to be expected 
for a variety of further experiments. 

We point out, that we are assuming adiabatic follow- 



ing of potential energy curves in all simulations. Non- 
adiabatic transitions are unlikely in the strong fields in- 
side the deflector, since the number of avoided crossings 
and their energy gaps generally increase with electric field 
strength. Moreover, the probability for non-adiabatic fol- 
lowing depends on the rate of change of the field strength, 
which is only due to the slow translational motion of 
the molecules. However, non-adiabatic transitions have 
been observed in different Stark decelerator beamlines 
at real [49] and avoided crossings [50] for small electric 
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FIG. 7: (Color online) Schematic illustration of the polariza- 
tion state of the YAG and the probe pulse with respect to 
the static electric field and the detector plane used to charac- 
terize alignment. The dashed line represents the propagation 
directions of the laser beams. Included is also a sketch of 
the resulting molecular alignment. Repeller, extractor and 
ground refers to the electrostatic plates of the velocity map 
imaging spectrometer. 



fields. Similarly, when the deflected molecules in the ex- 
periments reported here enter a fleld-free region, scram- 
bling of population over the various M components of 
their rotational state will occur. 

For small molecules, like OCS or CICN, the prepara- 
tion of an ensemble of molecules, all in a single quantum 
state will be feasible with the present setup. For these 
systems, the number of quantum states that are popu- 
lated in a supersonic jet is signiflcantly smaller compared 
to large asymmetric top molecules like IB or BN. The 
spacing between neighboring quantum states is larger 
and the number of avoided crossings smaller. Thus, the 
differences in the effective dipole moment between in- 
dividual quantum states are larger and, therefore, the 
degrees of their deflection will vary considerably. 



B. Laser-Induced Alignment of 
Quantum-State-Selected Molecules 

We now turn to studying alignment induced by the 
YAG pulse. The basic experimental observables are 2D 
1+ ion images recorded when the iodobenzene molecules 
are irradiated with both the YAG pulse and the probe 
pulse. The geometry of the laser pulse polarizations with 
respect to the velocity map imaging spectrometer (VMI) 
is illustrated in Figure 7. The YAG pulse is linearly po- 
larized along the vertical direction, i.e., in the detector 
plane. The probe pulse is linearly polarized perpendicu- 
lar to the detector plane, which ensures that there is no 
detection bias on the molecular orientation in that plane. 
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FIG. 8: (Color online) 1+ ion images illustrating align- 
ment, recorded when the probe pulse Coulomb explodes the 
iodobenzene molecules seeded in He. The polarizations of 
the YAG and the probe pulses are kept fixed as illustrated 
in Figure 7. The labels "no deflection", "deflection" and 
"depletion" correspond to images recorded at lens position 
= 0.0 mm, 1.0 mm and -0.9 mm respectively, the latter 
two with the deflector at 10 kV (see Figure 5). The inten- 
sities of the YAG and probe pulse are 8 x 10^^ W/cm^ and 
5 X 10^^ W/cm^, respectively. The color scale indicates the 
relative number of ions. This color scale is the same for all 
subsequent flgures showing ion images. 



This results in a circularly symmetric 1+ image, when 
only the probe pulse is used (Figure 8 Al). When the 
YAG pulse is included (Figure 8 A2-A4) the 1+ images 
exhibit strong angular confinement along the polarization 
of the YAG pulse. 

The 1+ ions appear as two pairs of radially localized re- 
gions, corresponding to two different fragmentation chan- 
nels of the Coulomb explosion. The radius of the out- 
ermost (and weakest) pair of rings is approximately 
times larger than the radius of the innermost (and bright- 
est) pair of rings. Since the radius is proportional to the 
velocity of the ions the 1+ ions from the outermost pair 
of rings originate from a Coulomb explosion channel that 
releases twice as much kinetic energy as the channel pro- 
ducing the 1+ ions in the innermost pair of rings. As 
pointed out in several previous studies from our group 
(see for instance reference 51) this is only consistent with 
the innermost pair of rings originating from iodobenzene 
being doubly ionized by the probe pulse and fragmenting 
into an I++ CqR^^ ion pair, and the outermost pair of 
rings originating from 1+ ions formed by triple ionization 
and fragmentation into an I++ CqB.^^^ ion pair. The 
pronounced angular confinement observed in images A2- 
A4 is quantified by calculating the expectation value of 
(cos^ ^2d), where O2D is the angle between the YAG pulse 
polarization and the projection of the 1+ recoil velocity 
vector onto the detector plane. In this paper, (cos^ O2D) 
values are calculated only from ions detected in radial 
region corresponding to the I++ CqK^^^ channel. By 
doing so, the YAG intensities probed are restricted to a 
narrow range close to the maximum value, as the high 
nonlinearity of the multiphoton process occurs efficiently 
only in the spatial regions close to the focal point of the 
YAG beam. 

Image Al, recorded with only the probe pulse present, 
should correspond to a target of randomly oriented 
molecules. As expected the image is circular symmet- 



8 



2.3 X 10^° W/cm^ 2.3 x 10^^ W/cm^ 1.2 x 10^^ W/cm^ 




FIG. 9: (Color online) I^ion images illustrating alignment 
at different intensities of the YAG pulse, recorded when the 
probe pulse Coulomb explodes iodobenzene molecules seeded 
in 20 bar Ne. The labels "no deflection" and "deflection" 
correspond to images recorded at lens position y = 0.0 mm 
(deflector turned off) and 2.15 mm (deflector at 10 kV), re- 
spectively. The intensity of the probe pulse is 5 x 10^^ W/cm^ . 



ric and (cos^^2d) = 0.515. When the YAG pulse is in- 
cluded (image A2) a pronounced angular confinement is 
observed along the polarization of the YAG pulse and 
(cos^ 6>2D)is increased to 0.947. These observations are in 
complete agreement with previous studies [38] . When the 
deflector is turned on and the laser foci moved to the edge 
of the most deflected molecules (at the position marked 
in Figure 5 a), corresponding to molecules in the lowest 
rotational states, the angular confinement is further en- 
hanced (image A3) leading to a (cos^ 02d) value of 0.968. 
By contrast, when the experiment is conducted on the 
least deflected molecules in the depleted region (image 
A4), corresponding to molecules in the highest rotational 
states, the alignment is weakened and (cos^ ^21)) =0.900. 

We repeated the alignment measurements when 
iodobenzene was seeded in Ne. The results are displayed 
in Figure 9. Like in the He case a pronounced improve- 
ment is observed when deflected rater than undeflected 
molecules are employed. Figure 9 shows images of 1+ 
recorded at three different intensities of the YAG laser 
for both undeflected and deflected molecules seeded in 
Ne. The effect of the deflector is clearly seen when com- 
paring, for instance, image Bl (deflected) and Al (unde- 
flected). At this low YAG intensity (2.3 x 10^^ W/cm^) 
weak alignment is obtained in the non-deflected beam 
with (cos^^2d)= 0.695. Going to the edge of the de- 
flected molecular beam (position indicated in Figure 5 
b ) a clear enhancement is observed (image Bl) with 
the (cos^6>2d) value rising to 0.869. Also, at high YAG 
intensity (1.2 x 10^^ W/cm^) the difference in angular 
confinement comparing the undeflected (image A3) to 
the deflected molecules (image B3) is visible. While 
(cos^ ^21)) =0.929 represents the limit of the degree of 
alignment of iodobenzene seeded in Ne in the undeflected 
beam, emloying the deflector leads to an unprecedented 
degree of laser-induced alignment of {co^^ 620) = 0.972. 
We note that {cos^ 620) = 1 would correspond to the 
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FIG. 10: Degree of alignment as a function of the YAG in- 
tensity for iodobenzene seeded in 20 bar Ne. The labels "no 
deflection" and " deflection" correspond to images recorded at 
lens position ^ = 0.0 mm (deflector turned off) and 2.15 mm 
(deflector at lOkV), respectively. The intensity of the probe 
pulse is 5 X 10^^ W/cm^ 



quantum mechanically unfeasible situation of perfectly 
ID aligned molecules. 

To quantify the angular information of the images, 
{cos'^ 02d) values are plotted as a function of YAG in- 
tensity, the results are displayed in Figure 10. Even at 
very low laser intensities a high degree of alignment can 
be obtained from an ensemble of quantum-state-selected 
molecules. The tendency shown in this graph with a 
steep rise and early saturation of the degree of alignment 
agrees with previous results investigating the dependence 
of alignment on the rotational temperature of the en- 
semble of molecules [38]. Effectively, the quantum-state 
selection corresponds to a "colder" albeit non-thermal 
beam (see Section III A). 

Note that the contrast between the undeflected and 
the deflected beam is expected to be greater if Ne is used 
instead of He as a carrier gas. In the undeflected beam 
the maximum degree of alignment that can be achieved 
is smaller in Ne because the rotational cooling in the 
supersonic expansion is less effective due to the lower 
stagnation pressure [37]. Additionally, in the deflected 
beam a better degree of alignment is expected for Ne, 
as the efficiency of the quantum-state selection in the 
present setup is significantly enhanced due to the longer 
residence time in the defiector (see Section III A). 



C. Laser-Induced Orientation of 
Quantum-State-Selected Molecules 

Next, we discuss orientation due to the combined ac- 
tion on the molecules by the YAG pulse and the static 
electric field (Egtat) from the VMI electrodes [5, 6]. 
Figure 11 illustrates the polarization state of the YAG 
and the probe pulse with respect to the static electric 



9 



repeller 



extractor ground 




YAG 



detector 



^^probe 



FIG. 11: (Color online) Schematic illustration of the polariza- 
tion state of the YAG and the probe pulse with respect to the 
static electric field and the detector plane used to character- 
ize orientation. The dashed line represents the propagation 
direction of the laser beams. Also, a sketch illustrating the 
molecular orientation is included. 



Al probe 

• 




A2 80° 

( 
1 

d 


A3 ^ 70° 

h 


A4 60° 

I 


A5 45° 

I 


A6 30° 

s 


Bl ^ 90° 

3 




B2 _ 100° 

f 


B3 ^ 110° 

? 

8 


B4 120° 

9 


B5 135° 

I 


B6 150° 

8 


CI probe 

1 




C2 80° 
t 


C3 ^ 70° 

m 

8 


C4 60° 
• 

d 


C5 45° 
« 

Q 


C6 30° 

6 


Dl 90° 
o 




D2 100° 
t 

t 


D3 110° 

s 

X 


D4 120° 

8 


D5 135° 




D6 150° 



FIG. 12: (Color online) I^ ion images illustrating orienta- 
tion for different values of /3, recorded when the circularly 
polarized probe pulse Coulomb explodes the iodobenzene 
molecules seeded in 90 bar He. The labels "no deflection" 
and "deflection" correspond to images recorded at lens posi- 
tion y = 0.0 mm (deflector turned off) and 1.0 mm (deflector 
at 10 kV), respectively. The intensity of the YAG and the 
probe pulse is 8 x 10^^ W/cm^ and 5 x 10^^ W/cm^, respec- 
tively. Estat = 594 V/cm. 



field of the VMI electrodes. The important difference 
compared to the alignment data is that the YAG polar- 
ization is rotated away from the axis perpendicular to 
the static field. Thus, the orientation data result from 
geometries where the angle /3, between the YAG polar- 
ization (the C-I bond axis) and the static electric field, is 
different from 90°- see Figure 11. To image orientation a 
circularly polarized probe pulse is used. This ensures that 
any molecule will be ionized - and thus detected - with 
the same probability independent of /3. This circularly 
polarized probe will induce some bias on the angular dis- 
tribution of the 1+ ions (see Figure 12 Al and CI), how- 
ever, importantly, it is up/down symmetric. Figure 12 
shows 1+ ion images for different /3 values for both de- 
flected and undeflected molecules seeded in He. As men- 
tioned the circularly polarized probe alone gives rise to 
an image that exhibits some angular confinement with 
(cos^ ^2i:))=0.70 (Figure 12 Al). Consequently, including 
the YAG pulse at /3 = 90° results in an image (Figure 12 
Bl) that appears slightly different from the correspond- 
ing image with a linearly polarized probe (Figure 8 A2), 
but still shows that the molecules are tightly aligned. 

Focussing first on the non-deflected data of Figure 12 
(row A and B) two prominent changes are observed as 
the polarization of the YAG pulse is gradually rotated 
away from the detector plane (images A2-A6 and B2- 
B6). First, the location of the 1+ rings shifts closer to 
the center of the images. This is due to the fact that the 
C-I axis alignment, and thus the emission direction of 
the 1+ ions, follows the YAG pulse polarization. When 
the C-I axis is aligned at an angle P the magnitude of 
the 1+ velocity vector recorded on the detector will be 
reduced by the factor sin(/3). The detrimental effect on 



the radial (velocity) resolution is obvious ai /3 = 135° /45° 
(image A5 and B5) and 30°/150° (image A6 and B6), 
where the two 1+ explosion channels, I++ C6H5+ and 
I++ C6H5^+, become indistinguishable as they merge in 
the 2D projection onto the detector plane. 

Secondly, as the YAG pulse polarization is turned away 
from 90°the up/down symmetry of the images, charac- 
teristic for the alignment data described in Section III B 
(and Figure 12 column 1), is broken. For images with 
90° < < 180° (images B2-B6) more 1+ ions are de- 
tected in the upper part, whereas for 0° < < 90° (im- 
ages A2-A6) more 1+ ions are detected in the lower part. 
The asymmetry becomes more pronounced as the YAG 
polarization is rotated closer to the axis of the static field. 
We interpret these observations as orientation due to the 
combined effect of the YAG laser field and the projec- 
tion of the static electric extraction field (Egtat) on the 
YAG polarization axis. This projection (numerical value: 
I cos(/3) I -E 

Stat increases as (3 is rotated towards 0° or 
180°, which is expected to cause an increase of the orien- 
tation [5, 6], in agreement with the experimental findings. 

As discussed in Section HI A all states of iodobenzene 
are high-field seeking, hence the orientation is expected to 
place the I-end of the molecules towards the repeller plate 
(see Figure 11), where the electrical potential is highest, 
because the dipole moment of iodobenzene is directed 
along the C-I axis pointing from iodine ( "negative end" ) 
towards the phenyl ring ("positive end"). The expected 
resulting molecular orientation at a given angle of /3 is 
shown in Figure 11. Thus, for 0° < < 90° the 1+ 
ions are expected to preferentially be ejected downwards, 
and for 90° < /3 < 180° they will be ejected upwards. 
This is in agreement with the up/down asymmetry in 
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FIG. 13: (Color online) I^ion images of deflected iodobenzene 
seeded in 20 bar Ne, illustrating orientation at two different 
intensities of the YAG pulse. The images are recorded at lens 
position y = 2.25 mm with the deflector at 10 kV. The inten- 
sity of the probe pulse is 5 x lO^'^ W/cm^. Estat = 594 V/cm. 



the images. 

The alignment (images Dl) and orientation (images 
C2-C6 and D2-D6) improves significantly when the de- 
flector is turned on and the foci of the lasers are moved 
to the position of the most deflected molecules [position 
marked in Figure 5 a)]. The markedly better orientation 
resulting in a much more pronounced up/down asymme- 
try is clearly visible even when the the YAG pulse is only 
turned slightly away from perpendicular, i.e., comparing 
deflected and undeflected images for /3=100° (image C2 
and A2) and (3 = 80° (image D2 and B2). From the 
previous discussion it appears that the highest degree of 
orientation is achieved when /3 is rotated towards 0° or 
180°. This is clearly seen from the images in row C and 
D and, again, the improvement obtained with deflected 
molecules is striking - compare image C6 to A6 (or D6 
to B6). 

Similar orientation measurements were conducted for 
iodobenzene seeded in Ne instead of in He. Figure 13 
shows 1+ images at a series of P values for two differ- 
ent intensities of the YAG pulse recorded with the de- 
flector at 10 kV at the position marked in Figure 5 b). 
Compared to the images displayed in Figure 12 row A 
and B, although recorded at slightly different intensities 
of the YAG pulse, a significant improvement is observed. 
Furthermore even at low intensity of the YAG a high 
degree of orientation for iodobenzene seeded in Ne is 
achieved. At the same time some loss in the angular 
confinement, i.e., in the alignment degree, is visible. We 
assign the clear improvement in the up / down asymmetry 
to the more stringent state selection in Ne compared to 
He as described in Section HI A, hence, manifesting itself 
in a higher degree of orientation. 

To quantify the up/down asymmetry, i.e., the de- 
gree of orientation, we determine for each image the 
number of 1+ ions, N(I+)up, in the upper part of the 
I^ + CeHs^ and 1+ + CeHs^^ channels (i.e., ions de- 




FIG. 14: Orientation of iodobenzene, seeded in either He or 
Ne, represented by the number of I^ ions in the upper half 
of the image (Nup) divided by the total number of I^ ions in 
the image (Ntotai) as a function of f3. For the experiments 
conducted with He as a carrier gas, curves C (lens position 
y = 1.0 mm and deflector at 10 kV) and D (lens position 
2/ = 0.0 mm and deflector off), the intensity of the YAG pulse 
was 7.8 X lO"*^"*^ W/cm^. For the Ne carrier gas, curves A 
and B (lens position y = 2.25 mm and deflector at 10 kV), 
the intensities of the YAG pulse are displayed in the insert. 
The intensity of the probe pulse for all curves was fixed at 
5 X 10^^ W/cm^ Estat = 594 V/cm. 



tected in the upper half of the images) as well as the to- 
tal number of ions, N(I )totai (= N(I+)up + N(I+) 

down I • 

This ratio, as a function of is displayed in Figure 14. 
Focussing first on curves C and D, representing iodoben- 
zene seeded in He, the difference between the data for the 
deflected molecules and the data obtained with the de- 
flector turned off is striking and shows the advantage of 
selecting the lowest-lying rotational states for strongly in- 
creasing the degree of orientation. The further improve- 
ment when Ne is used instead of He is clear from curves A 
and B. These two curves also show that the pronounced 
degree of orientation is maintained when the intensity 
of the YAG pulse is lowered by an order of magnitude 
compared to the maximum value of 1.2 x 10^^ W/cm^ 



IV. CONCLUSION AND OUTLOOK 

In conclusion, we have shown that deflection of cold 
molecular beams with an inhomogeneous static electric 
field enables the selection and the spatial separation of 
the most polar quantum states, i.e., the lowest-lying 
rotational states. The method demonstrated here is 
complementary to state-selection for small molecules 
using a hexapole focuser, which has been suggested 
to be applied for improved alignment and orientation 
experiments [52] and recently been experimentally 
demonstrated [53]. While a hexapole focuser only works 
for small molecules in low- field-seeking quantum states. 
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beam deflection will apply broadly to a wide range of 
molecules, from diatomics to large biomolecules. The de- 
flection is strongest for molecules with a large permanent 
dipole moment to mass ratio. For a given molecule the 
deflection is optimized by employing stronger deflection 
fields, increasing the length of the deflector, or lowering 
the speed of the molecule, for instance, by using neon 
rather than helium as a carrier gas. For small molecules, 
the preparation of an ensemble of molecules all in a 
single quantum state should be feasible. As an appli- 
cation of the state-selected molecules we showed that 
selection of iodobenzene in low lying rotational states 
allows to achieve unprecedented degrees of laser-induced 
adiabatic alignment and mixed laser- and static-field 
orientation. In particular, we demonstrated that strong 
alignment and orientation can be maintained even when 
the intensity of the alignment pulse is lowered to the 
10^^-10^^ W/cm^ range. This can reduce unwanted 
disturbance from the laser field in future applications 
of adiabatically aligned or oriented molecules. We note 
that it should be possible to improve the degree of 
orientation obtained here simply by increasing the static 
electric field. Due to experimental constraints this was 
not implemented in the present work. 

Getting access to cold molecules in the gas phase 
typically involves using a molecular beam from a su- 
personic expansion that usually consists of more than 
99 percent carrier gas and less than one percent of 
the specific molecules. In several types of experiments 
the atomic carrier gas can contribute to, or even 
completely overshadow the particular signal measured. 
The electrostatic defiection naturally separates the 
polar molecules from the unpolar carrier gas and thus 
removes this unwanted background. We conservatively 
estimate the density in the original molecular beam to 
be 10^^ molecules/cm^. In the experiments presented 
here the density in the defiected part of the beam is 
approximately 10^^ molecules/cm^. We foresee this 
density to be suflicient for a variety of applications, such 
as photoelectron spectroscopy with VUV, EUV [54], 
or x-ray light sources, including attosecond pulses, or 
high harmonic generation experiments with fs laser 
pulses [55]. For all of these applications the separation 
of the molecular target from the carrier gas might be of 
great relevance. 

Additionally, brute-force orientation [1, 2], that is, the 
spatial orientation of polar molecules using strong dc 
electric flelds, will beneflt from the state-selected sam- 
ples similar to what was demonstrated here. In fact, 
the states that are deflected the most are also oriented 
the most in a dc electric fleld. To illustrate the achiev- 
able orientation, we have calculated the ensemble av- 
eraged orientation in a homogeneous electric fleld of 
250 kV/cm for iodobenzene for a thermal ensemble of 
1 K. For this ensemble {cosO) = 0.757 is obtained. For 
a deflected, quantum-state-selected sample of iodoben- 



zene molecules at 1 % of the undeflected peak intensity 
(see Table I), an increased ensemble averaged orienta- 
tion of {cosO) = 0.905 is predicted. For nonadiabatic 
alignment [4, 56] it was demonstrated, that the dynam- 
ics and, importantly, the degree of alignment and orien- 
tation depends strongly on the initial rotational state dis- 
tribution [57-61]. Selection of rotational states is, there- 
fore, highly advantageous for nonadiabatic laser-induced 
schemes to control the spatial orientation of molecules, 
as recently demonstrated for NO molecules. [53] In gen- 
eral, complete elimination of the rotational tumbling of 
an asymmetric top molecule requires, that all three prin- 
cipal axes are conflned along laboratory flxed axes. This 
is the area of 3D alignment and orientation, subjects that 
have been treated theoretically and experimentally [62- 
67]. A straightforward extension of the current work 
would be signiflcantly improved 3D alignment and ori- 
entation of polar molecules by employing state-selected 
molecules and we have recently performed such experi- 
ments [67]. 

For large (bio-)molecules, typically multiple structural 
isomers (conformers) [68] are present even at the low tem- 
peratures in a supersonic jet [69]. These conformers of- 
ten exhibit large, and largely different dipole moments 
what can be exploited to spatially separate them using 
inhomogeneous electric flelds. This separation has re- 
cently been demonstrated in an alternating gradient fo- 
cusing (AG) selector for the cis- and trans-conformers 
of 3-aminophenol [7]. However, also the static fleld of a 
deflector can be used to spatially isolate individual con- 
formers [70]. In general, the ability to achieve very high 
degrees of alignment and orientation is of great inter- 
est for a number of applications. We therefore believe 
state-selected molecules could be very beneflcial for areas 
such as photoelectron angular distributions from flxed-in- 
space molecules [71], (ultrafast) diffraction with electron 
or x-ray sources [8, 9] and time-resolved studies of light- 
induced stereochemistry [72]. Furthermore, Janssen and 
coworkers have shown that the ability to select a single 
rotational state with a hexapole focuser enables new pos- 
sibilities for studying directional dynamics of fragments 
in photodissociation of small molecules [26]. The deflec- 
tion method strongly increases the number of molecules 
to which single rotational state selection, and subsequent 
orientation, can be applied. In particular, it will offer 
access to studies of photoinitiated processes in oriented 
targets of larger asymmetric tops. 



APPENDIX A: SIMULATION OF BEAM 
PROFILES 

In this Appendix, the simulations of the experimen- 
tal deflection proflles are described in detail. First, 
the calculation of the Stark effect for asymmetric top 
molecules is explained in Section Al. Section A 2 de- 
scribes how single-quantum-state deflection proflles are 
obtained from trajectory simulations using the calculated 
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Stark curves. Finally, the single-quantum-state profiles 
are averaged in a suitable way to simulate a thermal en- 
semble for a given rotational temperature, which is de- 
scribed in Section A3, and fitted to the experimental 
data (Section A 4). 

1. Stark Effect Calculations 

In order to calculate the adiabatic energy curves for 
asymmetric top molecules, the Hamiltonian matrix is set 
up in the basis of symmetric top wavefunctions. In the 
presence of an electric field, only M is a good quantum 
number for an asymmetric rotor. J, which is a good 
quantum number in the field-free case, is mixed by the 
field, whereas K is mixed by the molecular asymmetry. 



{JKM I i^rot I JKM) 
{JK ± 2M I i^rot I JKM) 



{JKM I fia I JKM) 
(J + IKM I I JKM) 



{JK ± IM I /i5 I JKM) 
( J + li^ ± IM I /i5 I JKM) 



{JK ±1M \ fic\ JKM) 
{J^1K±1M \ JKM) 



For the correct assignment of the states to the "adiabatic 
quantum number labels" ^ M, i.e., to the adiabati- 
cally corresponding field-free rotor states, one has to clas- 
sify the states according to their character in the electric 
field symmetry group [75, 76]. This symmetry classifica- 
tion can be performed by applying a Wang transforma- 
tion [77] to the Hamiltonian matrix. If the molecule's 
dipole moment is along one of the principal axes of iner- 
tia, the matrix will be block diagonalized by this trans- 
formation according to the remaining symmetry in the 



Thus, the adiabatic energy curves can be calculated for 
the different M levels individually by setting up and di- 
agonalizing the M matrices including all J and K lev- 
els. An accurate description of higher rotational quan- 
tum states also requires including centrifugal distortion 
constants. The Hamiltonian H of an asymmetric rotor 
molecule with dipole moment j~L in an electric field of 
strength E can be written as the sum of the Hamiltonian 
i^rot of an asymmetric rotor in free space in Watson's 
A-reduction [73] and the contribution due to the Stark 
effect i^stark as 

H = i^rot + ^Stark (Al) 

Following references 73 and 74, the corresponding matrix 
elements are: [81] 



(A2) 



(A3) 
(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 
I 

field and the blocks can be treated independently. For 
arbitrary orientation of the dipole moment in the iner- 
tial frame of the molecule the full matrix must be diago- 
nalized. In any case, this process ensures that all states 
(eigenvalues and eigenvectors) obtained from a single ma- 
trix diagonalization do have the same symmetry, and, 
therefore, no real crossings between these states can oc- 
cur. Therefore, by sorting the resulting levels by energy 
and assigning quantum number labels in the same order 
as for the field-free states of the same symmetry yields 



B + C 



{J{J+1)-K^)+AK^ 



Ajj\j + If - AjkJ{J + l)K^ - AkK" 



B-C 



5jJ{J + l) 



^/J{J + 1) - K{K ± 1)^/J{J + 1)-{K± 1){K ± 2) 
MK 



J(J+1) 
{JKM \iia\J + IKM) 



v/(J+l)2-i^V(J + l)2-M2 

(J + l)v/(2J + l)(2J + 3) 
M^{JtK){J±K + 1) 
2J{J+l) 
{JK±1M I fib I J+IKM) 



, V(J±i^ + l)(J±X + 2)V(J+l)2-M2 
2(J+l)V(2J+l)(2J + 3) 



My'{JTK){J±K+l) 

2J{jTT) '^^^ 

{JK±1M I fic I J+IKM) 
.^{J±K + 1)(J ±K + 2)^(J + 1)2 - M2 
2(J + l)v/(2J + l)(2J + 3) 
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the correct adiabatic labels. These calculations are per- 
formed for a number of electric field strengths - typically 
in steps of 1 kV/cm from kV/cm to 200 kV/cm - and 
the resulting energies, I^stark(£^), are stored for later use 
in simulations using the libcoldmol program package [78] . 

2. Monte Carlo Simulations of a Molecular Beam 

Simulations of the electrostatic beam deflection in our 
setup are performed with the home-built software pack- 
age libcoldmol [78]. Trajectories for individual molecules 
in a given rotational quantum states are obtained from 
numerical integration of the 3D equations of motion using 
a Runge-Kutta algorithm. The initial phase space distri- 
bution of the molecular packet in the transverse spatial 
coordinates, x and is described by the mean values and 
widths of circular uniform distributions. For the veloc- 
ity coordinates, Gaussian distributions characterized by 
their mean values and full widths at half maximum are 
used. Also the initial time spread of the molecular beam, 
which corresponds to the opening time of the valve, is 
described by a Gaussian distribution centered around to- 
From the initial phase space distribution, the position of 
a molecule in phase space is randomly chosen. Then, it 
is propagated through the beamline, which includes all 
mechanical apertures of the experimental setup. From 
the electric field in the electrostatic deflector, which is 
calculated in two dimensions using finite element meth- 
ods (Comsol Multiphysics 3.4), and the Stark energy (see 
Appendix Al), the force acting on the molecule in the 
transverse directions is obtained: F = — VWstark(^)- 
The electric field is taken to be constant along z. Finally, 
the deflected molecules are propagated through field-free 
space to the detection region. For each quantum state, 
10^ trajectories are calculated yielding single-quantum- 
state deflection profiles Is{y)- 

3. Deflection Proflles 

The spatial deflection profile for an ensemble of 
molecules at a given rotational temperature, /(y. Trot), is 
calculated from the single-quantum-state deflection pro- 
files, /s(y), as follows: 

1 ^ 

/(^,Trot) = - V^.(Trot)/.(^) (AlO) 

w ^-^ 

s = l 

Here, N is the number of quantum states included in the 
simulations and Ws{Trot) is the population weight for a 
given quantum state: 

Ws(T,ot) = 9M9nse ^^rot (All) 

with Wo being the field-free potential energy of the 
ground state and Ws the field-free energy of the cur- 
rent state; = 1 for M = and gM = 2 otherwise; 



^ns accounts for the nuclear spin statistical weight of the 
current state. For both, benzonitrile and iodobenzene, 
^ns = 5 for Ka even and ^ns = 3 otherwise. The normal- 
ization is given hy w = XlfLi '^s- 

In some cases, a molecular beam cannot be accurately 
described by a single rotational temperature [46] . There- 
fore, a two-temperature model is implemented in our 
simulations. In the two-temperature model, the low- 
temperature part of the molecular beam is calculated 
as described above. In order to realistically simulate 
the high-temperature component in the molecular beam, 
a very large number of quantum states would have to 
be included in the simulations for large asymmetric top 
molecules. Because these Monte Carlo simulations are 
very time consuming, the high-temperature component 
is approximated by adding a small fraction of undeflected 
molecules to the deflected beam profile. This approxima- 
tion is well justified, because typical rotational tempera- 
tures for the high-temperature component in a molecular 
beam are on the order of 10 K and, for these temper- 
atures, most of the molecules reside in high rotational 
quantum states that have a small Stark shift and, there- 
fore, remain almost undeflected. Thus, the deflection 
profile, /2T(y, ^rot, <?), in the two-temperature model is 
calculated as 

hAy. Trot, q) = q- i{y. Trot) + (1 - ^) • iud{y) (A12) 

where < q < 1 and Iud{y) denotes the undeflected spa- 
tial beam profile that is obtained when both electrodes 
are grounded. 



4. Fit of Rotational Temperature 

In order to calculate the rotational temperature of the 
molecular beam, the simulated deflection profiles are fit- 
ted to the experimental data using a Nelder-Mead sim- 
plex algorithm. In the fitting procedure, the difference 
between simulated and measured deflection profiles is 
minimized, where the rotational temperature Trot of the 
low-temperature component, the fraction q of molecules 
in this low-temperature component, and a general inten- 
sity scaling factor s of the deflected profiles (with respect 
to the undeflected beam profile) are used as fitting pa- 
rameters. All deflection curves measured for different 
high voltages are fitted simultaneously to determine the 
rotational temperature. 
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